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We study the resistivity vs. electric field dependence p{E) of a 2D hole system in SiGe close to the 
B = metal-insulator transition. Using p as a "thermometer" to obtain the effective temperature 
of the holes Te{E), we find that the p{E) dependence can be attributed to hole heating. The 
hole-phonon coupling involves weakly screened piezoelectric and deformation potentials compatible 
with previous measurements. The damping of the Shubnikov-de Haas oscillations gives the same Te 
values. Thus the p{E) dependence and the i5-field "scaling" do not provide additional evidence for 
a quantum phase transition (QPT). We discuss how to study, in general, true i5-field scaling and 
extract the ratio of the QPT characteristic lengths. 



PACS numbers: 71.30.+h, 63.20.Kr, 73.43.Nq 

The recent observations of a metal-insulator transition 
(MIT) in two-dimensional (2D) electron or hole systems 
in high mobility silicon metal-oxide-semiconductor field 
effect transistors (Si-MOSFETs) ||,|l and in certain het- 
erostructures P-B[ have triggered important experimen- 
tal and theoretical efforts to understand the unexpected 
metallic behavior ^^ ■ This consists in a decrease of the 
resistivity p for decreasing temperature T, at densities Ps 
larger than the critical density Pc- In contrast, dp/dT < 
on the insulating side of the MIT {ps < Pc)- The ex- 
trapolated finite resistivity at T=0 is in conflict with the 
scaling theory of localization for non-interacting particles 
which predicts localized states in 2D systems llfl]. The 
MIT occurs for values of r-g (the ratio of a carrier pair 
Coulomb energy to the Fermi energy) much larger than 
one, suggesting the existence of a non-Fermi liquid due 
to strong interactions. The MIT could thus be the signa- 
ture of a T=0 quantum phase transition (QPT) between 
new metallic and insulating phases. However, for some 
systems, e.g. p-SiGe, it seems that the positive dp/dT is 
due to quasi-classical processes masking the usual weak 
localization which should appear again at low enough 
temperatures P,[ll|-p^ . At present, the nature of the 
MIT is still a subject of ongoing discussions. 

A striking feature of the MIT is that at low temper- 
ature, the dependence of p as a function of the electric 
field E is similar to p{T), i.e. dp/dE < for ps < Pc 
and dp/dE > for ps > Pc [|0,|l4|-|l3l ■ The physical ori- 
gin of this observation is an open question. Irrespective 
of the possible microscopic explanations of the p{E,ps) 
dependence pjlia], a QPT implies that p scales with E 
if it scales with T [091: close to the critical point, two 
characteristic length scales are associated with T and E, 
L^{T) - T-i/^ and Le{E) ~ ^-i/(^-Hi)^ and p depends 
only on the ratio of the smallest of them to the correla- 
tion length 1^ '^ 1(5„| {z is the dynamical exponent, v 



the correlation length exponent, and (5„ — (ps —pc)/Pc)- 
Thus p{T,ps) (for ^ ^ 0) depends only on {Sr^/T^''-" 
on each side of the transition, and p{E,ps) (at low tem- 
perature) depends only on \5n\ /E^^'^^'^^^^ . The scaling 
analysis of p{T,ps) and p{E,ps) allows a separate ex- 
traction of z and v. E scaling has been observed in sev- 
eral systems exhibiting T scaling [|7|,p^-p7[ . However, a 
natural question is whether the p{E) dependence results 
from carrier heating |n9[. At low temperature, the weak 
carrier-phonon coupling can lead to an effective temper- 
ature of the carriers, Te{E), larger than the lattice tem- 
perature Ti |p|,p0|-p8[. Thus, the p{E) dependence can be 
due to the Te{E) dependence of p. 

In this paper, we demonstrate experimentally that for 
a 2D hole system in p-SiGe exhibiting the MIT features 
with T and E scaling, the resistivity vs. £'-field de- 
pendence can be attributed to hole heating. We find 
p{E,ps) = p{T = Te{E),Ps) with a law T^{E) compati- 
ble both in shape and magnitude with known electron- 
phonon coupling models and data. The same T^i^E) is 
obtained from the damping of the Shubnikov-de Haas 
(SdH) oscillations. Thus, in the present case, the p{E,ps) 
dependencies and the _E-field "scaling" are not arguments 
in favor of the QPT interpretation of the "MIT" and the 
separate extraction of v and z is not possible. As heating 
effects depend on the semiconductor structure, while E 
scaling is a crucial issue for the QPT, we discuss the con- 
ditions required to study it for various systems in spite 
of carrier heating. We show how the ratio L^/Le can be 
obtained experimentally. 

The experiments were performed on Sio.85Geo.15 quan- 
tum wells sandwiched between undoped Si layers [f7|,p^. 
The 2D hole gas was formed in the triangular potential 
well at the Si/SiGe interface located on the boron doped 
side. It occupies the lowest heavy-hole subband, with an 
effective mass of about 0.25mo. Two gated Hall bars (SI 
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FIG. 1. (a)-(b) Resistivity as a function of temperature 
(for E — > 0), at different densities for sample SI. (c)-(d) 
Resistivity as a function of electric field for various densities 
and lattice temperatures. 



and S2) were used. Their width W is 100 /ini, and their 
length L (between voltage probes) 233 jjlth for SI, and 
125 /im for S2. By varying the gate voltage, their den- 
sities can be tuned between 0.50 and 1.75 x 10^^ cm"^. 
The mobilities at 200 mK increase with ps from 600 to 
5500 cmVVs (resp. 1000 to 7400 cmVVs) for SI (resp. 
S2). An ungated sample, with a density of 3.9 x 10^^ 
cm^^ and a mobility of 7800 cm^/Vs has also been used. 
The temperature range was 70 mK - 1.4 K. To check that 
the temperature T given by the thermometer fixed in the 
copper sample holder was that of the lattice Ti , the low 
current resistance of another Hall bar etched onto the 
same substrate was used as a thermometer. Its temper- 
ature remained close to T whatever the current in our 
samples, p — {V / I){W / L) and E — V/L were obtained 
from the current / and the voltage drop V between the 
voltage probes, using a four point DC technique, with a 
current (10 - 300 pA for low E measurements) periodi- 
cally reversed at a frequency 0.03 - 0.3 Hz. 

The MIT features appear in the p{T,ps) dependence 
for ^ ^ [Fig. 0(a)-(b)]. dp/dT < is found for 
Ps < Pc~l-3 X 10" cm-2 (r, « 6), while dp/dT > 
for Ps > Pc, at temperatures above a threshold which 
decreases when ps increases. At the largest densities 
dp/dT is positive. The p{T,ps) data plotted as a func- 
tion of \5n\/T'' for T > 0.3 - 0.5 K collapse on two 
branches (not shown), demonstrating a scaling behav- 
ior with b = 0.45 ± 0.04. The values 0.35 and 0.62 
[g|,^ were found in p-SiGe. The MIT characteristics and 
i?- field scaling appear also in the p{E,ps) curves (Fig. g), 
Pc remaining unchanged. When p is plotted as a func- 
tion of \Sn\ /E" (a=0.19 ± 0.02), the curves faU on two 
branches (inset of Fig. 0). They contain only a part of 
the data (symbols on Fig. S) : removing the low E points 
is justified by Le < i$ (assuming a real QPT), while for 
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FIG. 2. Resistivity of sample SI as a function of electric 
field at various hole densities: from ps — 0.92 x 10^^ cm~^ 
(top curve), to ps — 1.75 x 10^^ cm~^ (bottom curve). The 
lattice temperature is 120 mK. The lines give the whole set 
of data for _E > 10 mV/cm, while the symbols correspond to 
the fraction of the data used in the curves presented in the 
inset. Inset: resistivity as a function of |(5„| /i?"'^^. 



large E, it can be related to a microscopic limit on Le- 
Figure |l| shows that the general shape and the minima 
of p{E) are the same as those of p(T), suggesting hot 
hole effects. Their contribution is studied as follows: (i) 
the p{E) dependence is assumed to result from the hole 
temperature rise from T/ to T^. due to Joule heating; (ii) 
for each value oi E, p is used as a "thermometer" giving 
Tg as the temperature T at which the same value of p 
is measured for E —^ 0; (iii) the relationship between Tg 
and the power per carrier Pe = VI /{psWL) = E'^ /{pps) 
is used to study the validity of the hot carrier assumption 
(i). Figure H(a) shows the PE{Te) dependence for various 
lattice temperatures and for densities around Pc- 

The power loss between a degenerate 2D electron sys- 
tem and the lattice 3D acoustic phonons is given by pG] 



Pe = A{T^ - Tf ) + A'{Tf - Tf '). 



(1) 



The first term corresponds to the deformation potential 
coupling, with a — 5 (resp. 7) for weak (resp. strong) 
screening. The second term is piezoelectric coupling, to 
be considered for our coherently strained SiGe samples 
pl| , p9[ , with a' — 3 (resp. 5) for weak (resp. strong) 
screening. The exponents may be decreased by one in 
disordered systems because of "dynamic" rather than 
"static" screening |23]. A and A' are prefactors related 
to the deformation (S„) and piezoelectric (cpz) coupling 
constants [|o|,|l). To test the vahdity of Eq. 1^) for 
our measurements, we write it: Pe + {ATf + A'T" ) = 



{AT^ + A'Tf). Figure |(b) shows that by adding to 
Pe a constant Po{Ti) chosen separately for each T;, the 
whole set of curves for a given density falls on the same 
master curve. Hence PE{Te,Ti) + Po{Ti) depends only 
on Te, and this dependence is the sum of two power 
laws. We have verified that the Po{Ti) dependence is 
the same, thus proving that our data are in agreement 
with Eq. (m. Similar results are obtained for all densi- 
ties and samples. For ps close to pc, a power law with 
a' « 3 (resp. a « 5 — 6) dominates at low (resp. 
large) Tf,. These power laws are confirmed by fitting 
the PsiTe) data with Eq. (g). Fig. |(a) shows the good 
quality of the fit when a — 5 and a' = 3 are imposed. 
For Ps > 1.35 X 10" cm'^, a term A"(T^ - Tf) is 
added in the fit, corresponding to the cooling through 
the contacts which increases when p decreases [Q. A" 
has the same order of magnitude as the value given by 
the Wiedeman-Franz law. Our data are thus compati- 
ble with Eq. (H), the two main cooling processes being 
hole-phonon coupling through weakly screened deforma- 
tion and piezoelectric potentials. Weak screening at low 
temperatures has been pointed out in p-SiGe [ pi]j2^ ] and 
Si-MOSFETs ||,||. A and A' are obtained from a fit 
with Eq. (|]) assuming a — b and a' — i [Fig. ||(a)]. Fol- 
lowing Ref. Q we then obtain S^ = 2.7 ± 0.3 eV and 
ep2 = (3.4 ± 0.9) X 10-3 C/m^ for all the densities. Our 
Sti value is compatible with the S^ « 3.0 eV found in 
p-^iGe at ps = (3.5 - 13) x lO^^ cm^^ ^M. Our Cp^ 
is somewhat lower than the Cpz ~ 1.6 x 10 C/m^ ob- 
tained in p-Sio.8Geo.2 [0, however in Rcf. |2^] Cpz was 
found to be much smaller than in Ref. pij . 

In order to further prove the validity of the heating 
analysis, the Te extracted at i? = is compared to the 
temperature obtained using the damping of the SdH os- 
cillations |2|,|2|,||,|||. To get rid of the p^^iT) depen- 
dence due to carrier scattering |2^ , the "thermometer" is 
the difference between a minimum oi p^xiB) and the pre- 
vious maximum. Its T dependence results from a differ- 
ent physical situation (the density of states oscillations) 
than at B = 0. The magnetic field values (_B w 1 T) 
are low to ensure that the electron-phonon coupling laws 
extracted are close to those at i? = |^j|. As shown 
in the inset of Fig. ^(a) for ps — 1.19 x 10^^ cm"^ and 
Ti = 150 niK, the two methods are in very good agree- 
ment. A similar agreement is obtained for other Ti values 
and forps = 3.9 x 10^"'^ cm~^. Thus, the p{E) dependence 
can be attributed to hole heating, implying that it does 
not bring new physical information on the possible MIT 
when compared to the p{T) dependence. 

To study £'-field scaling in spite of heating effects, the 
condition Le{E) < Lt],[Te{E)] must be fulfilled. Figure 
shows how i$ and Le depend on E. The vertical scale 
is arbitrary because L^ and Le are not given by the 
QPT theory, but assuming z — I for strongly interacting 
particles [|^, L* - T^^ and Le ^ E-^l"^. Although 
the "metallic" p{T) in our samples is compatible with 
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FIG. 3. (a) Power per hole injected in the 2D system 
as a function of the holes effective temperature, for two 
gated samples and three densities, at several lattice tem- 
peratures (arrows). For clarity, the curves corresponding 
to Ps = 1.19 X 10" cm-2 (resp. to p, = 0.55 x 10" 
cm"^) are shifted upwards by 100 (resp. by 10*). The 
lines corresponding to the two lower densities are fits to 
Pe = A{Tl - rf ) -I- A!{T't - Tf)\ while for the larger density 
a third term yl"(r^ — Tf) has been added in the fit formula. 
Inset: Pe as a function of the effective temperature extracted 
using the p thermometer at _B = (o) or the damping of 
the SdH oscillations (-I-). (b) The sum of the power per hole 
and a constant -Po(31) chosen for each Tj, as a function of the 
hole effective temperature, for various lattice temperatures 
Ti = 100 mK (O), 150 mK (A), 200 mK (v), and 380 mK 
(D). The lines give the slopes of the power laws T^ and T^. 



a Fermi liquid description |12| , it is instructive to con- 
sider the possibility of a QPT: the full line corresponds 
to L^\Tq(E)\ in our case, for ps k. p^,, using the experi- 
mental law Te{E), obtained in the E interval indicated 
by the symbols. An upper limit for E arises since the 
longitudinal potential drop V = EL has to be kept small 
compared to the gate voltage. As we attributed the mea- 
sured p{E) dependence to heating, the Le{E) line (long 
dashed) crosses the L^\Te{E)\ curve beyond this interval. 
The E-fie\A scaling can be investigated only beyond this 
crossing point, ("i? scaling" interval in Fig. 0), provided 
L$ and Le do not reach their microscopic limit. 
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FIG. 4. The logarithm of L* and Le as a function of the 
electric field. The full line represents L,f,[Te{E)] calculated us- 
ing Eq. (§) at ps = 1.31 x lO" cm"^ and Ti = 100 mK. The 
symbols correspond to our experimental points. Short-dashed 
line: L.,,[T^{E)] estimated for Si-MOSFETs close to the MIT. 
Long-dashed line: Le{E) assumed to be the same for SiGe 
heterostructures and Si-MOSFETs. The vertical lines sepa- 
rate the hot carrier region from the T and E scaling regions 
for p-SiGe. The grey area gives the limits above which the 
_B-field scaling is improved in Ref. [17] 

The short dashed line in Fig. | is L^[Te{E)] for Si- 
MOSFETs, estimated using recent measurements of the 
power loss [|4[|o). Again, the Lb(^) < Lq.[Te{E)] 
prescription leads to a lower limit for E. This agrees 
with the results of Ref. ||l^ where the quality of the in- 
field scaling is improved when E is larger than a min- 
imum value (the grey area in Fig. El corresponds to 
their cut-offs of 50 to 500 pW). For p-GaAs, the a = 5 
exponent quoted in Ref. p8| leads to a similar situa- 
tion, but the thermal coupling should be larger than in 
Si-MOSFETs. Unscreened piezoelectric coupling would 
lead to an upper £'-ficld limit. A better knowledge of 
the power loss laws would allow a quantitative use of the 
Le{E) < L^{Te) prescription. The field Ec defined by 
Le{Ec) — L^\Te{Ec)\ could be extracted as the limit 
beyond which the experimental PE(Te) law differs from 
the power loss law, thus yielding the important physical 
result L^[T)/Le{E) = [r,(£;,)/r]V-(^/^^)i/(-+i). 

The ratio of the exponents in the experimental scal- 
ing laws p(|5n| /T^) and p(|5n| /E"") has been proposed 
as an indicator of the p{E) dependence origin |1^ . For 
the QPT scaling, a/h — z/{z -I- 1) = 0.5. For carrier 
heating, Eq. (0) gives T^ ~ (£;2/")p-i/"^ thus a/b = 2/ a 
neglecting A' , Ti and the p{T) dependence. We find ex- 
perimentally a/b = 0.42 ± 0.08, while 0.4 and 0.67 are 
expected for our a and a', thus the a/b criterion can 
hardly be used ||3l| ]. 

In summary, we have shown that in a 2D hole system 
exhibiting the B — MIT characteristics, with T and 
E scaling, the p{E) dependence close to the MIT could 



be interpreted as being due to hole heating. Thus, in 
our case, the experimental p{E,ps) and £'-field scaling 
are not an indication that the MIT features can be at- 
tributed to a QPT, and the separate extraction of v and 
z is not possible. However, there is an E interval de- 
fined by Le{E) < L^[Te{E)] where i?- field scaling can 
be investigated independent of hot carrier effects. 

The authors thank T. Jolicocur for helpful discussions. 
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We used this equation with Us 
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For weak screening the conductivity dependence of the 
prefactor given in Ref. [E3] should not be relevant. 
Similarly, in Si-MOSFETs, unscreened electron-phonon 
coupUng [|^J25| would lead to 2/a = 0.4. 



